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ABSTRACT 

We report on 8.7 and 7.6 yr of Rossi X-ray Timing Explorer (RXTE) obser- 
vations of the Anomalous X-ray Pulsars (AXPs) RXS J170849.0-400910 and 
IE 1841—045, respectively. These observations, part of a larger RXTE AXP 
monitoring program, have allowed us to study the long-term timing, pulsed flux, 
and pulse profile evolution of these objects. We report on four new glitches, one 
from RXS J170849.0-400910 and three from IE 1841-045. One of the glitches 
from IE 1841—045 is among the largest ever seen in a neutron star in terms of 
fractional frequency increase. With nearly all known persistent AXPs now seen 
to glitch, such behavior is clearly generic to this source class. We show that in 
terms of fractional frequency change, AXPs are among the most actively glitch- 
ing neutron stars, with glitch amplitudes in general larger than in radio pulsars. 
However, in terms of absolute glitch amplitude, AXP glitches are unremarkable. 
We show that the largest AXP glitches observed thus far have recoveries that are 
unusual among those of radio pulsar glitches, with the combination of recovery 
time scale and fraction yielding changes in spin-down rates following the glitch 
similar to, or larger than, the long-term average. We also observed a large long- 
term fractional increase in the magnitude of the spin-down rate of IE 1841—045, 
following its largest glitch, with Az//z> = 0.1. These observations are challeng- 
ing to interpret in standard glitch models, as is the frequent occurence of large 
glitches given AXPs' high measured temperatures. We speculate that the stellar 
core may be involved in the largest AXP glitches. Furthermore, we show that 
AXP glitches appear to fall in two classes: radiatively loud and radiatively quiet. 
The latter, of which the glitches of RXS J170849.0-400910 and IE 1841-045 
are examples, show little evidence for an accompanying radiative event such as a 
sudden flux increase or pulse profile change. We also show, however, that pulse 
profile and pulsed flux changes are common in these AXPs, but do not apprear 
closely correlated with any timing behavior. 
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Introduction 



The past decade has seen significant progress in our knowledge of the o bservational 
properties of Anomalous X-ray Pulsars (AXPs; see IWoods fc Thompson! 120061 ; iKaspil 120071 . 
for recent reviews) . F rom a timing point of v i ew, the presence of binary companions has been 
practically ruled out ( Mereghetti et al. 1998 : Wilson et al. 1999b . and sub sequently their po- 
tential for great rotational stability was demonstrated (IKaspi et al.lll999l). thereby allowing 



the discovery that AXPs c an exhibit spin-up glitches (IKaspi et al. 



2000 



2003 



DairOsso et al.l 120031 1 . and large (factor of ~10) torque variations (jGavriil fc Kaspi 



Kaspi fc Gavriil 



2004 ). From a radiative point of view, AXPs are now known to show a var iety of differ- 



ent variabil ity phenomena, in c 



like bursts (Gavriil et al 



Ibrahim et al. 



2004 



2002 



Israel et al 



2004 



Woods et a 



2007a 



Dib et al 



uding long-lived flares (jGavriil fc Kaspil l2004r). short SGR 



2005) , large outbursts (IKaspi et al.l 12003 



ted), and slow, low-level flux and pulse profile variability (IDib et al 



2007b: Tarn et al.l 2006: Tarn et al, submit- 



2007a: Gonzalez et al. 



submitted). Spectrally , though previously studied only in the soft X-ray band, A XPs are now 
seen in the radio band ( Camilo et al. 2006h. through the mid- (Wang et al. 2006h and near-IR 



Israel et al. 



Rea et al 



2002 



2004: 



2002 



Wang fc Chakrabartv 



2002; Hulleman et al 



2004 



Tarn et al 



2004 



Durant fc van Kerkwijkl 120051 ) , in the optical range (e.g., iKern fc Martin 



Dhillon et al.l 120051 ) , up to hard X-ray energies (IKuiper et al.l 120061 ) . The evidence 



thus far argues strongly that AXPs, like their close cousins, the Soft Gamma Repeaters, are 
magnetars - young, isolated neutron stars powered by a large magnetic energy reservo ir, 
with surface fields of > 10 14 — 10 15 G (jThompson &: Duncan! Il996l : iThompson et al.ll2002l ). 



In spite of this progress, however, many aspects of AXPs remain mysterious. Particularly 
so are their variability properties. What is the origin of the variety of different types of 
variability? Although bursts can be exp lained as sudden crustal yields, slower evolution (e.g., 
Gavriil fc Kaspil 2004 ; IDib et al. 2007a ) has been suggested to be due to slow magnetospheric 
twists (IThompson et al.ll2002l ). Some support for this pic ture has been argued to come from 



observed correlations between fl ux and spectral hardn ess (IWoods et al.ll2004l : iRea et al.ll2005 



Campana et al.ll2007l ). although lOzel fc Guverl (120071 ) argue that such a correlation need not 
originate uniquely from the magnetosphere and could be purely thermal. At least some 
radiative variability has been seen to be correlated with timing behavior. The best example 
of this occured in the 2002 outburst of AXP IE 2259+586 in which the pulsar suffered a 
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large spin-up glitch app a rently simultaneous ly with a major X-ray outburst fjKaspi et al. 
20031 ; IWoods et al.ll2004l ). Ilsrael et al.l (j2007al ) describe a similar radiative outbur st in AXP 



CXOU J164710. 2— 455216, and report a large contemporaneous glitch, as did iDib et al. 
J2007bh recently for AXP IE 1048.1-5937. By contrast, AXP RXS J1708 49.0-400910 



exhib i ted two glitches with no eviden ce for a corresp o nding radiative event (IKaspi et al. 



2000l ; IKaspi fc Gavriill 120031 ) although iDafTOsso et al.l (|2003f) suggested po ssible low-level 
pulse profile changes associated with the second glitch. ICampana et al.l (120071 ) also suggested 
that observed flux and spectral changes may be associated with glitches and predicted a 
third glitch would be observed after mid-2005 on the basis of an observed flux increase and 
apparently correlated spectral changes. 

Here we report on 8.7 and 7.6 yr of monitoring of RXS J170849.0-400910 and IE 1841-045, 
respectively, using the Proportional Counter Array (PCA) aboard the Rossi X-ray Timing 
Explorer (RXTE). We report the discovery of one new glitch and three new glitch candi- 
dates in RXS J170849. 0-400910 as well as three new glitches in IE 1841-045, including 
one of the largest glitches, in terms of fractional frequency increase, thus far observed in 
any neutron star. We also present pulsed flux time series for RXS J170849. 0—400910 and 
IE 1841—045 which reveal little or no evidence for correlated changes with glitches, al- 
though RXS J170849. 0—400910 shows low-level pulsed flux variability at many epochs. We 
also report a pulse profile evolution analysis which shows that both pulsars' profiles are 
evolving slowly with time, though in neither case does this evolution show a clear correla- 
tion with timing behavior. These results demonstrate that AXPs RXS J170849. 0—400910 
and IE 1841—045 are frequent glitchers. They also demonstrate that although AXP timing 
glitches can occur simultaneously with significant long-lived radiative enhancements, they 
need not always do so. 



2. Observations 

The results presented here were obtained using the PCA on board RXTE. The PCA 
consists of an array of five collimated xenon/methane multi-anode proportional counter 
units (PCUs) operating in the 2—60 keV rang e, with a total effec tive area of approximately 
6500 cm 2 and a field of view of ~1° FWHM ( Jahoda et al.lll996l ). Our 294 observations of 



RXS J170849. 0-400910 and our 136 observations of IE 1841-045 are of various lengths (see 
Tables [T] and [2]) • Most were obtained over a period of several years as part of a long-term 
monitoring program, but some are isolated observations (see Figures [1] and [2]) . 

For the monitoring, we used the GoodXenonwithPropane data mode except during Cy- 
cles 10 and 11 when we used the GoodXenon mode. Both data modes record photon arrival 
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times with l-/is resolution and bin photon energies into one of 256 channels. To maximize 
the signal-to-noise ratio, we analysed only those events from the top xenon layer of each 
PCU. 



3. Phase-Coherent Timing 

Photon arrival times at each epoch were adjusted to the solar system barycenter. Result- 
ing arrival times were binned with 31.25-ms time resolution. In the RXS J170849. 0—400910 
timing analysis, we included only events in the energy range 2—6 keV, to maximize the 
signal-to-noise ratio of the pulse. Similarly, for IE 1841—045 we included events in the 
energy range 2—11 keV. 

Each barycentric binned time series was folded using an ephemeris determined iteratively 
by maintaining phase coherence as we describe below. Resulting pulse profiles, with 64 phase 
bins, were cross-correlated in the Fourier domain with a high signal-to-noise template created 
by adding phase-aligned profiles from all observations. The cross-correlation returned an 
average pulse time of arrival (TOA) for each observation corresponding to a fixed pulse 
phase. The pulse phase <p at any time t can usually be expressed as a Taylor expansion, 

(Pit) = Mto) + v {t - to) + ^o(t - to? + \^it - t ) 3 + • • (1) 

2 6 

where v = 1/P is the pulse frequency, v = du/dt, etc., and subscript "0" denotes a parameter 
evaluated at the reference epoch t = to- The TOAs were fitted to the above polynomial using 
the pulsar timing software package TEMPO0. 

Note that we also searched for X- ray bursts in e ach 2- 20 keV barycentered, binned time 



series using the methods described in iGavriil et al.l (120041 ). however no bursts were found in 



any of our RXS J170849.0-400910 or IE 1841-045 data sets. 



3.1. Timing Results for RXS J170849.0-400910 

Figure [3] and Table [3] summarize our results for RXS J170849. 0—400910. The pulsar's 
spin evolution can be characterized by steady spin-down, punctuated by sudden episodes of 
spin-up, i.e., glitches, in addition to candidate glitch events and apparently random noise. 
We provide in Table Opulse ephemerides for inter-glitch ranges labelled as in the top panel of 



http: / / www.atnf.csiro.au/ research/pulsar /tempo 
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Figure [3j Residuals after subtraction of these models are shown in the next panel of Figure [3j 
Overall the models describe the data well. However, particularly when our timing precision 
was highest (i.e., before 2003), some low-level but significant deviations are seen on time 
scales of weeks to months. Their o rigin is unknown but is likely related t o "timing noise," 



commonly seen in oth er AXPs (e.g-.lKaspi et a. 



in radio p ulsars (e.g . jArzoumanian et al. 



1994 



1999 



Gotthelf et a 



Hobbs et al. 



2004 



200 21) and ubiquitousl y 



Livingstone et al.ll2005l ). 



Note that iDib et al.l (l2007al ) performed simulations which showed that pulse profile changes 
similar to those observed in this source (see §4.ip did not result in timing offsets significantly 
larger than our reported TOA uncertainties. Hence, the features in the timing residuals 
reported here are not a result of pulse profile changes. 

In addition to the two previously reported glitches (which we have reanalysed, finding 
resul t s consistent with thos e already in the literature; see lKaspi et al.ll2000l ; iKaspi fc Gavriil 
20031 ; iDairOsso et al.ll2003l ). we have identified a third unambiguous glitch that occured near 
MJD 53551 (2005 June 30). Note that the exact glitch epoch is unknown due to our non- 
continuous monitoring; we report an epoch for which the phase jump is zero. This is because 
a non-zero phase jump at the time of the glitch would suggest an unphysically large torque on 
the star. This third glitch had fractional frequency jump Avjv = 2.7 x 10 -6 , and no obvious 
recovery. This glitch amplitude is intermediate between those of the previous two observed 
glitches, and the lack of recovery is similar to what was seen in the first glitch, but in marked 
contrast with the second glitch, as is clear from Figure HI A sudden change in post-glitch 
spin-down rate for the third glitch is difficult to constrain, because of a possible additional 
glitch that occured not long after, as we describe below. Indeed glitc h-induced long-term 
changes in i> aside from that following the first glitch, as described by IKaspi et al.l (120001 ). 
are difficult to identify given the apparent timing noise processes. Table H] summarizes the 
parameters of the three certain glitches of RXS J170849. 0—400910, assuming a glitch model 
consisting of a permanent change in v and i> and a frequency change v& that decayed on a 
time scale of r^, i.e., 



v 



u Q (t) + Au + Au d e 



+ AV (t - tg 



(2) 



where u (t) is the frequency evolution pre-glitch, Av is a instantaneous frequency jump, Av d 
is the post-glitch frequency increase that decays exponentially on a time scale T d , t g is the 
glitch epoch, and Ai> is the post-glitch change in the long-term frequency derivative. 

For the second glitch, residuals after subtraction of a simple glitch with fractional ex- 
ponential recovery have clear remaining trends, as is clear in Figures [3] (second panel) and H] 
(b ottom panel). System atic trends after simple glitch model subtraction were also reported 
bv IWoods etaD ((200J) for the 2002 glitch in IE 2259+ 586. We also find this in the largest 
glitch in IE 1841-045 (see §E£2). IWoods et all J2004T ) showed that for IE 2259+586, the 
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glitch fit was significantly improved by adding an exponential growth term. We have tried 
fitting this model to the second glitch from RXS J170849. 0—400910 but find no improvement, 
with the preferred growth term consistent with zero. 

In addition to the new certain glitch we report above, we find strong evidence for an 
additional three glitches, each having fractional amplitude similar to the first certain glitch 
seen in this source. The properties of these candidate glitches are summarized in Table [5j 
Timing residuals around the epochs of these glitches are shown in Figure EJ in the top panel. 
Residuals following the subtraction of a glitch model are shown in the middle panel of that 
Figure. We refer to these as candidates only because a 4-th order polynomial fit to the 
same data results in similar residuals (bottom panel of Fig. |5j Table [6]) without the need 
to invoke a sudden event. The distinction between true gl itches and timing noise is often 



difficult to make for small-amplitude glitches, as discussed by lKaspi et al.l (120001 ) . One way to 
distinguish, at least statistically, is that apparent discontinuities attributable to timing noise 
should not have a preferential direction, i.e., apparent spin-down 'glitches' should be seen 
too. An examination of the frequency panel in Figure [3] reveals apparent frequency jumps at 
the candidate glitch levels in both directions, suggesting one or more of the candidates could 
indeed be timing noise. Continued monitoring to acquire a larger database of such apparent 
discontinuities will help clarify this issue. 



Subsequent to our submission and posting of this paper, llsrael et al.l (l2007bl ) posted the 
results of a similar analysis of a subset of these same data. Some of their results are consistent 
with ours however others differ. They reported two glitches, the first of which corresponds 
to our second candidate glitch (Table EJ). For that glitch, the reported fit parameters are 
similar though not identical to ours. Their second glitch corresponds to our third glitch in 
Table HI For that glitch, the reported frequency jump at the glitch epoch was similar to ours 
but the jump in frequency derivative was significantly different. We find that this difference 
is due to their inclusion of more post-glitch TOAs when fitting the glitch. We did not include 
these TOAs because of a candidate event that occurs shortly thereafter, but which Israel 
et al. did not report. In addition to this difference, the frequency value reported in their 
post-glitch ephemeris, 0.09088624(2) Hz, is 42a away from the value 0.090885327(8) Hz that 
we measure at the same epoch using our post-glitch ephemeris. The numbers in parentheses 
are lex uncertainties. We do not understand this difference. 



3.2. Timing Results for IE 1841-045 



Figure [6] and Table [7] summarize the long-term timing behavior of IE 1841—045. As for 
RXS J170849. 0-400910, the spin evolution of IE 1841-045 is well characterized by regular 
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spin-down punctuated by occasional sudden spin-up events, plus timing noise. Ephemerides 
in Table [7] are given for the glitch-free intervals indicated in the top panel of Figure [6) As for 
RXS J170849. 0—400910, the long-term timing residuals show some unmodelled trends wh ose 



origin is unknown. We consider these trends timing noise, as did iGotthelf et al.l (120021 ) in 
their analysis of ~2 yr of data from this object. Note that the ephemeris in Table [7] labeled 
B2 is the same as that labeled B except for the omission of data immediately post-glitch (see 
caption to Fig. [6]). 

The frequency panel in Figure [6] and first panel in Figure[7Jmake clear that IE 1841—045 
suffered a large glitch, with significant recovery, near MJD 52460 (2002 July 5). This epoch 
is estimated, as for all glitches reported in this paper, by taking the epoch at which the phase 
jump is zero. Note that for this glitch there were several such epochs and the one we are 
reporting gives the most conservative frequency jump assuming an exponential recovery. The 
least conservative possible frequency jump is ~ 50% larger. The glitch fractional amplitude 
was Avjv = 1.6 x 10~ 5 (see Table [8]), among the largest yet seen from any neutron star. A 
fraction Q = Av^j (Az/^ + Au) = 0.64 of the glitch recovered on a time scale of 43 days. This 
glitch is thus similar to the second certain glitch seen in RXS J170849. 0—400910, and to the 
2002 glitch in IE 2259+586, which also showed significant recoveries on time scales of weeks. 
Also, like the second glitch of RXS J170849.0-400910, this large glitch in IE 1841-045 is 
not well modelled by Equation [2], as is clear in the residuals plot in Figure [7J Accompanying 
this frequency glitch was a substantial long-term increase in the the magnitude of i>, with 
fractional increase Ai> jv = 0.0959 ± 0.0007. This is discussed further in §6.11 

Because of the sparsity of the data around the glitch epoch, we found an alternate 
ephemeris for the period of time covered by ephemeris B (Fig. [6j Table [7]). The fit parameters 
are v = 0.0849041677(17) Hz, v = -2.852(3) x 10" 13 Hz s _1 , v = -2.47(8) x 10~ 21 Hz s~ 2 , 
and d 3 v/dt 3 = 8.8(7) x 10~ 29 Hz s~ 3 at the reported glitch epoch MJD 52464.00448, with 
RMS phase residual of 0.019. This ephemeris disagrees with ephemeris B in the shape of 
the recovery (see dotted curve in Panel 3 of Fig. [6]) but agrees with it after the end of the 
recovery. Using the parameters of this alternate ephemeris, the change in v at the glitch 
epoch would be 2.20(3) x 10~ 7 Hz much smaller than the one reported in Table [HJ However, 
we hesitate to interpret the glitch using this ephemeris because of the very unusual and 
unique shape of the recovery it predicts. Note that this alternate ephemeris also shows a 
long-term increase in the magnitude of i> after the glitch. 

We also report the detection of two additional, smaller glitches, as summarized in Table [8] 
and displayed in Figures [6] and [7J Neither glitch displays significant recovery, and both are 
well modelled by a simple permanent frequency jump. 
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4. Pulse Profile Changes 

Another interesting AXP property we can study thanks to RXTE monitoring is the 
evolution of the pulse profile. We performed a pulse profile analysis on each AXP using 
FTOOLS version 5.3. 1q We used the following steps: for each observation, we ran the 
FTOOL make_se to combine the GoodXenon files. We then used the FTOOL f asebin 
to make a phase-resolved spectrum of the entire observation with 64 phase bins across the 
profile. When we ran f asebin, we selected layer 1 of the detector, disregarded the propane 
photons, and included the photons from PCUS 1, 2, 3, and 4. We omitted PCU 0, for 
which an indepen dent analysis of AXP 4U 0142+61 revealed spectral modeling irregularities 



(IDib et al.ll2007al ). f asebin also took care of barycentering the data. For each observation, 
we then used seextrct to make a phase-averaged spectrum for the same set of detector 
layers and PCUs. The phase-averaged spectrum was then used by the perl script pcarsp to 
make a response matrix. 

We loaded the phase-resolved spectra and the response matrices into the X-ray Spectral 
Fitting Package (XSPEC@) and selected photons belonging to three energy bands: 2—10, 
2—4, and 4—10 keV. Using XSPEC, we extracted a count-rate pulse profile for each of the 
energy bands. The profiles included XSPEC-obtained la error bars on each of the phase 
bins. To obtain a pulse profile in units of count rate per PCU, we divided the overall profile 
by a PCU coverage factor that took into account the amount of time each PCU was on. 

We then aligned the 64-bin profiles with a high signal-to-noise template using a similar 
cross-correlation procedure to the one used in the timing analysis. Then, for each glitch-free 
interval, we summed the aligned profiles, subtracted the DC component, and scaled the 
resulting profile so that the value of the highest bin is unity and the lowest point is zero. 



4.1. Profile Analysis results of RXS J170849.0-400910 

Average profiles for RXS J170849. 0—400910 in the three energy bands are presented in 
Figure [HJ In a given band, the different profile qualities are due to different net exposure 
times. Energy dependence is clearly visible to the eye as well as small fluctuations. For 
example, in the 2-4 keV band, the small peak off the main pulse has clearly fluctuating 
intensity. This small peak gets larger at higher energy, as seen in the 2-10 keV band. In the 
4-10 keV band, the smaller peak seems to blend with the main low-energy peak to yield a 



2 http: / /heasarc. gsfc.nasa.gov/ftools 
3 http://xspec. gsfc.nasa.gov Version: 11.3.1 
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broad single peak structure, although fluctuations in that structure are apparent. 

To study these fluctuations quantitatively, we subjected each profile to a Fourier anal- 
ysis. Figure M shows the evolution of the first three profile harmonics with time. Although 
there are hints of variation in all energy bands, only variations in the hard 4-10 keV band are 
statistically significant (as determined by the x 2 statistic from a fit to a constant value); the 
decline of the second and third harmonics in the hard band have probabilities of 0.0007% 
and 0.0012%, respectively, of being due to chance. Thus in the hard band the profile is 
certainly becoming more sinusoidal, in agreement with what is inferred by eye. 

The above analysis shows that the profile is changing, but not whether these changes 
are truly correlated with the glitch epochs, since changes could be occuring t hroughout. To 



search for pulse profile changes correlated with glitch epochs, as were claimed by lDall'Osso et al 



(120031 ). we divided glitch-free intervals into several sub- intervals (typically of duration ~30 
days) for which independent profiles were created. The number of sub-intervals was chosen 
by trading off signal-to-noise ratio for time resolution. These sub-interval profiles were then 
Fourier analysed. The evolution of the Fourier powers in the first three harmonics in the 
2-10 keV profile are shown in the top panel of Figure (TDJ 

To determine whether the apparent fluctuations are statistically significant, we fit a 
constant value to each time series and from the x 2 of the best fit, found that the probabilities 
of the fluctuations being due to random noise are 68, 96 and 69% for n = 1, 2, 3, respectively. 
This analysis thus shows no evidence for profile changes associated with the glitch epochs, 
including the second glitch. However the reduced signal-to-noise ratios in the sub-interval 
average profiles, required for interesting time resolution, makes us insensitive to subtle profile 
changes. To search for glitch-correlated pulse profile changes in a different way, for each sub- 
interval we calculated the reduced x 2 °f the difference between that sub-interval's average 
profile and the previous one. The time series of these x 2 values is shown in the bottom panel 
of Figure [TOl There is clearly no evidence for any profile change at the second glitch, or at the 
third certain glitch. There is some hint of profile changes at the first and second candidate 
glitches, however a K-S test shows that our x 2 values as a group have a probability of 39% of 
originating from x 2 distribution. Interestingly though, the probability of the single high x 2 
value we measure at the second candidate glitch occuring randomly is only 1.0 x 10~ 6 ; that 
at the first glitch is 1.7% and at the first glitch candidate is 0.4%. Thus we do find possible 
evidence in this analysis for glitch-correlated pulse profile changes, though the best evidence 
for significant changes occurs only at two candidate glitches, i.e., the lowest amplitude events. 
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4.2. Profile Analysis results of IE 1841-045 

Summed profiles for IE 1841—045 in three energy bands for the five glitch-free intervals 
defined in the top panel of Figure [6] are shown in Figure [TT1 As for RXS J170849. 0—400910, 
some low-level profile fluctuations are suggested, particularly in the relative amplitude of the 
leading and trailing sides of the large single peak in the 2-10 keV band (though clearly this 
peak could also be considered the blend of two or more adjacent peaks). 

As for RXS J170849.0-400910, we quantify the profile fluctuations of IE 1841-045 via 
Fourier analysis. Figure [12] shows the evolution of the first three profile harmonics with time 
in each energy band. Interestingly, in contrast to RXS J170849. 0—400910, here the profile 
changes are most prominent in the soft 2-4 keV band, in which the fraction of power in the 
fundamental of the profile in interval A2 (see Fig. [6]) decreased, then slowly relaxed back to 
the previous range. However, a x 2 test shows the probability of this behavior being due to 
random noise is 18%, too large to exclude. 

To look for pulse profile changes correlated with glitches, again, sub-intervals within 
glitch-free intervals were chosen and summed profiles computed and Fourier analysed. The 
evolution of the first three harmonics is shown in the top panel of Figure [13j Fluctuations 
are apparent although none is particularly remarkable at any of the glitch epochs, including 
the first and largest, and the time series for n = 1,2,3 are all consistent with a constant 
value. This argues again against correlated profile and timing anomalies in this source thus 
far. As a confirmation, as for RXS J170849. 0—400910, a difference profile was calculated for 
each sub-interval by subtracting that interval's profile from the preceding one. The \ 2 values 
of these difference profiles are shown in the bottom of Figure [131 no significant features are 
present. 

5. Pulsed Flux Time Series 

RXTE monitoring also allows the study of the evolution of the pulsed flux of these 
sources. To obtain a pulsed flux time series for RXS J170849. 0-400910 and IE 1841-045, we 
did the following. First, for each observation, we used a procedure similar to that described 
in §4] to make a count rate per PCU pulse profile (with 64 phase bins across the profile and 
excluding PCU 0) in the energy range 2-10 keV. The profiles included XSPEC-determined 
la error bars on the flux value in each of the phase bins. 
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The pulsed flux for each of the profiles was calculated using the following RMS formula: 



\ 



2V(K 2 + 6/)-K 2 + ^ 2 )), (3) 



k=l 

where is the k th even Fourier component defined as = jf^2f =1 Pi cos (2irki/N), a ak 2 is the 
uncertainty of ajt, bk is the odd k th Fourier component defined as bk = jf^2iLiPi sin (2nki/N), 
<Jb k 2 is the uncertainty of bk, % refers to the phase bin, N is the total number of phase bins, pi 
is the count rate in the i th phase bin of the pulse profile, and n is the maximum number of 
Fourier harmonics to be taken into account. We used n = 6 for both RXS J170849. 0—400910 
and IE 1841-045. 

Our method for estimating the pulsed flux F is equivalent to the simple RMS formula 

F = a/ Y^iLi {Pi ~P) 2 (where pi is the count rate in the 2 th phase bin of the pulse profile 
and p is the average count rate), except that we have subtracted the variances (to eliminate 
the upward statistical bias) and only included the statistically significant Fourier components. 
For a detailed discussion on pulsed flux estimates, see Archibald, Dib & Kaspi (in prep.). 



5.1. Pulsed Flux Time Series for RXS J170849.0-400910 



Our pulsed flux time series for RXS J170849. 0—400910 is shown in Figure [3] (bottom 
panel) and again in Figure HH Each data point represents the average of pulsed fluxes 
measured over ~1 month. There appear to be frequent low-level pulse flux variations in this 
source. Although our error estimates on the pulsed fluxes include only statistical uncertainty 
(i.e., we have made no effort to estimate systematic uncertainties), we are given confidence 
that the fluctuations seen e.g., near MJD 52000 are real, given how stable the pulsed flux of 
IE 1841—045 is in the same time interval (see §5.21) . 

There are no large increases in pulsed flux followi ng any of the glitches, unlike what wa s 
seen following the 2002 glitch of AXP IE 2259+586 kaspi et al.lboO.j Iwoods et al1l2004l ). 
However, there is a possible pulsed flux enhancement prior to the second glitch, and a dip 
following it. Given this and the lack of clearly associated pulse profile changes coincident 
with glitch epochs (see gl), the glitches of RXS J170849.0-400910 appear to be "quiet," 
in the sense that they seem unaccompanied by significant pulsed radiative change. This is 
discussed further in SHU 



Figure [14] shows phase- averaged fluxes in the 0.5-10 keV band as me asured using a 
variety of focussing X-ray telescopes (IRea et al.ll2005l ; ICampana et al.l 120071 ). Interestingly, 
while the reported phase- averaged flux varies considerably (by a factor of ~1.6), and in 
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concert with the photon index as measured in the conventionally used blackbody/power-law 
spectral model, the 2-10 keV pulsed flux remains relatively constant. This suggests that the 
pulsed fraction of RXS J170849. 0—400910 is precisely anti-correlated with total flux, in such 
a way as to keep the pulsed flux near constant. This is discussed further in §HJ 

The origin of the apparent low-level pulsed flux variations is not clear, given the apparent 
lack of correlation with the phase-averaged flux. As shown by Archibald et al. (in prep.), 
a changing pulse profile can affect an RMS-based pulsed flux estimator such as that in 
Equation [3J To verify that our measured pulsed fluxes were not influenced by the changing 
pulse profile of RXS J170849. 0—400910 (see §@|, we also found the pulsed flux using an 
estimator based on the area under the profile (after baseline subtraction), which is, by 
definition, insensitive to pulse profile changes. With this method, we obtained qualitatively 
similar results for the pulsed fluxes. 

5.2. Pulsed Flux Time Series for IE 1841-045 

Our pulsed flux time series for IE 1841—045 is shown in the bottom panel of Figure [6j 
Each data point represents the average of pulsed fluxes measured over ~1 month. Note the 
increased scatter after MJD 52700 is due to decreased effective integration time, a result of 
the reduction in the average number of operational PCUs. The measured pulsed fluxes are 
consistent with being constant, with their probability of being due to random fluctuations 
52%. There is no evidence for any pulsed flux change at the glitch epochs. Thus the glitches 
of IE 1841—045 appear to be "quiet," at least in pulsed flux, on time scales comparable to 
or longer than our sampling time. 

6. Discussion 

6.1. AXP Glitches 

We have now observed a sufficiently large sample of AXP glitches that we can make 
meaningful phenomenological comparisons with glitches in radio pulsars, a much better 
studied phenomenon. Detection of systematic differences in AXP and radio pulsar glitch 
properties would be interesting as it could signal structural differences between magnetars 
and conventional radio pulsars. 

Figure [151 shows the fractional and non-fractional amplitude distributions of radio pulsar 
and AXP glitches. As is clear from the figure, although the fractional glitch amplitudes of 
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AXPs are generally large by radio pulsar standards, the AXP absolute glitch amplitudes, 
more directly related to the angular momentum transfer during the glitch, are neither es- 
pecially large nor especially small. Thus, glitching in neutron st ars is clearly not correlated 
with frequency as some studies of radio pulsars have suggested (ILyne et al.ll2000l ). 



Given the spectacular radiative outburst contemporaneous with the large 2002 IE 2259+586 
glitch, we can speculate that larger angular momentum transfers that occur in radio pulsars 
could result in e ven more dramatic outbursts in affected AXPs, possibly like th ose seen in 



XTE 1810-197 fllbrahim et al.ll2004h and in the AXP candidate AX J1845-0258 flTam et al. 



2006 ). Indeed a recent X-ray burst observed from CXOU J164710.2-455216 flMuno et al. 
20071 ) has been claimed to be accompanied by a very large (Av/v ~ 6 x 10 -5 , Av ~ 6 x 10 -6 ) 



glitch (llsrael et al.l l2007af) . and AXP IE 1048.1—5937 recently exhibited a large glitch and 
flux increase (IDib et al.ll2007bl . Dib et al., in prep.). However the lack of any observed ra- 
diative change in IE 1841—045 around the time of its first observed glitch, which was over a 
factor of two larger than that in IE 2259+586 in terms of absolute frequency jump, argues 
against this idea. Clearly, the data are indicating that AXP glitches, even large ones, can 
be either radiatively loud or quiet. 



Glitch activity has been defined as 



1 x - Au 
At^~ : 



(4) 



where At is t he total observing span and the sum is over all glitches, and includes decaying 
components (IMcKenna &: Lyndll990l ). We refer to a g as fractional activity, since it involves 
the sum of fractional frequency changes. One can also define an absolute glitch activity, 



A, 



At ^ 



(5) 



where the sum is over the absolute frequency changes (e.g- JWang et al.ll2000l ). The quantities 



a g and A g , introduced for the study of radio pulsars, are approximately interchangeable for 
those objects, given that the range of frequencies encompassed by glitching radio pulsars 
is relatively small. By contrast, when considering AXPs and their much smaller rotation 
frequencies, a comparis on with radio pulsars for a g and A g are very different (see, e.g., 
Heyl fc Hernquistlll999l ). Also, for establishing the average amount of spin- up imparted to 
the crust over time, the total frequency increase at each epoch is relevant. However, in some 
instances, the quantity of interest is the unrelaxed portion of the glitch, i.e., the permanent 
frequency jump only. In general, for radio pulsars, Q is small so this distinction is not 
important. However for AXPs, given the paucity of glitches we have observed thus far as 
well as the fact that several, particularly the largest, of these have had large values of Q 
(e.g., Q — 1 for the second glitch seen in RXS J170849. 0—400910), the distinction between 
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including the total frequency jump and only the unrelaxed frequency jump is important. We 
choose here to remain with convention and include the total frequency jump at each epoch 
when calculating a g and A g , although this choice should be kept in mind. 

With 8.7 yr of monitoring of RXS J170849. 0—400910, we can now reasonably calculate 
glitch activity parameters for this source using Equations H] and [5j If only counting the 
unambiguous glitches, a g = 2.9 x 10~ 14 s -1 and A g = 2.5 x 10 -15 s~ 2 . Including candidate 
glitches only increases these numbers by ~20%. With three glitches in 7.6 yr, IE 1841—045 is 
evidently a very active glitcher as well. Its glitch activity parameters are a g = 7.9 x 10~ 14 s _1 
and A g = 6.7 x 10~ 15 s -2 . Indeed A g for IE 1841—045 is the highest glitch activity seen 
thus far in any neutron star, radio pulsar or AXP, to our knowledge. We also calculated a 
tentative glitch activity for AX P IE 2259+586, for which we have ob served two glitches, the 
well documented one in 2002 fjKaspi et al. 2003 ; Woods et al. 2004 ) and a second, smaller 
glitch that occured very recently and had fractional amplitude 8.5 x 10~ 7 and no recovery 
(Dib et al., in prep.). Using these events and given that we have observed this source with 
RXTE for 9.4 yr, we find a g = 1.7 x 10~ 14 s" 1 and A g = 2.4 x 10~ 15 s _1 . 

We can plot these activities as a function of pulsar age (as estimated via spin-down 
age vjlv) and z>; see Figure [161 Previous authors have noted inte resting correlations on 



these plots for radio pulsars (e.g., iLyne et al.ll2000i ; IWang et al.ll2000i ); these are seen in our 



plots as well. Note that upper limits for some ra dio pulsars of relevant ages fall well below 
the apparent correlations (e.g.. IWang et al.ll2000l ); we choose not to plot those because, as 



discussed in that reference, a single glitch of average size would bring them roughly in line 
with the correlation. Note that the radio pulsar outlier at small age and high z> in all plots 
is the Crab pulsar, long-known to exhibit few and small glitches. We also looked for a trend 
in a plot of a g or A g versus surface dipolar field (as estimated via 3.2 x 10 19 y/pp G) but 
found none. 

As a group, the AXPs do not especially distinguish themselves when either activity, a g 
or A g is plotted versus spin-down age, though they do increase the scatter. This suggests a 
universal correlation with spin-down age. The same is true of A g plotted versus v. However, 
interestingly, the AXPs as a group all stand out on the diagram of a g versus z> (Fig. [T6|) . 
such that for similar spin-down rates, their fractional glitch activities are much larger than 
in radio pulsars. 



Link et al.l (119991 ) argued that a g provides a strict lower limit on the fraction of the mo- 



ment of inertia of the neutron star that resides in the angular momentum reservoir (generally 
assumed to be the crustal superfluid) tapped during spin- up glitches, I res . They showed that 
Ires/Ic > va g/\v\ = G, where I c is the moment of inertia of the crust and all components 
strongly coupled to it, and G is a "coupling parameter." For radio pulsars, they argued 
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for a universal G that implies I res /I c > 0.014. It is interesting to ask whether this same 
apparently universal relationship holds for AXPs. Fig ure [T71 shows G plotted versus age for 
radio pulsars and AXPs. As is clear, the iLink et al.l (119991 ) relation seems to hold for the 
radio pulsars, even with increased glitch statistics. Also, AXPs RXS J170849. 0—400910 and 
IE 1841—045 lie among the radio pulsars, suggesting similar reservoir fractions. However 
the outlier point, IE 2259+586, has G = 0.25, much larger than the others. Admittedly, 
for this AXP, a g is estimated from two glitches only, with the 2002 event greatly dominat- 
ing, so the values are tentative. Still, the large G, if real, suggests that at l east ~25% o f 



the s tellar moment of inertia is in t he angular momen tum reservoir (see also IWoods et al. 



2004J). We note that the analysis of ILink et al.l (119991 ) ignores recovery, important for the 
2002 IE 2259+586 glitch, which dominates its a g . However in the 2002 glitch, the recovery 
fraction was only ~19%, so even accounting for recovery, G for IE 2259+586 is surprisingly 
high. 



As described by Kaspi et al.l fj2Q03h and IWoods et all (J200J), the 2002 IE 2259+586 
glitch was unusual when compared with those of radio pulsars. Specifically the combination 
of the recovery time scale and the large recovery fraction Q conspired to make the pulsar 
spin down, for over two weeks post-glitch, at over twice its long-term average spin-down 
rate. Although spin-down rate enhancements post-glitch are often seen in radio pulsars (e.g., 
Flanaganlll990l ). they usually amount to only a few percent. A remarkably large post-glitch 
spin-down rate enhancement was seen also in the second glitch of RXS J170849. 0—400910 
and the first observed glitch of IE 1841—045, though to a lesser degree than in IE 2259+586. 
Of course a much larger increase in spin-down rate post-glitch in RXS J170849. 0—400910 or 
IE 1841—045 could have been missed due to our sparse sampling. 

One way to quantify the enhanced spin-down more precisely is using Equation [2] at 
t = 0, and noticing that the instantaneous spin-down rate at the glitch epoch due to the 
exponential recovery is given by Au^/r. Comparing this quantity for the AXP glitches 
that show recovery with the pre-glitch time-averaged spin-down rate z>, we find that for 
IE 2259+586 Av d jr = (8.2±0.6)z>, Au d /r = (0.64±0.6)z> for RXS J170849.0-400910, and 
Av d jr = (0.75 ± 0.08)z> for IE 1841—045, all very large by radio pulsar standards. 

The increase in spin-down rate post-glitch, at least for radio pulsars, is generally at- 
tributed to a decoupling of a small percentage of the moment of inert ia of the star, usually 
presumed to be part of the crustal superfluid (e.g.. |Pines fc Alparlll985l ). with constant exter- 
nal torque. If the observed AXP recoveries and temporarily enhanced spin-down rates were 
interpreted in the same way, it would imply that very large fractions (ranging from 0.4 to 0.9) 
of the moment of inertia of the star decoupled at the glitch, much larger than the crustal su- 
perfluid is reasonably expected to comprise, for any interior equation of state. To avoid this 
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problem, IWoods et al.1 (12004 ) suggested that a pre-glitch rotational lag between the crus t 
and superfluid might have temporarily reversed at the glitch (see also lAlpar et al.1 120001 ) . 
Then the observed larger spin-down rate post-glitch would be due to the crust transferring 
angular momentum back to the superfluid in order to reestablish equilibrium. 

In glitches, the equilibrium angular velocity lag between the crust and more rapidly 
rotating crustal superfluid is thought to be the origin of glitches. This lag is proposed 
to de velop because the crustal superfluid's angular momentum vortices, in many models 



^e.g., lAlpar et al.l Il984al ). become pinned to crustal nuclei and hence are hindered from 
moving outward as the star's crust and associa ted components are slowed by the external 
torque. How this lag could reverse is puzzling. IWoods et al.l (120041 ) suggested that a twist 
of magnitude 10~ 2 rad of a circular patch of crust offset in azimuth from the rotation axis 
could result in sufficient spin-down of the crustal superfluid to account for the properties of 
the 2002 glitch in IE 2259+586. They noted f urther that such a twis t also produces X-rays 
of the luminosity observed in that outburst ([Thompson et al.1 120021 ). The absence of any 
significant radiative changes at the time of largest glitches in RXS J170849. 0—400910 and 
IE 1841—045 is thus problematic for the crustal twist, and hence lag reversal, model. We 
note that it has been argued independently that a similar suggested lag reve rsal between 
crust and crustal superfluid in the Vela pulsar is unphysical ( Jahan-Miril 120051 ) . 



We also note that the large and long-term increase in the magnitude of v following the 
large glitch in IE 1841—045 (see §3.21) is also interesting. lAlpar et al.l (119931 ) showed that, 
ignoring transient terms, I res /I c > Ai>/i>. For the large IE 1841—045 glitch, this implies 
I res /I c > 0.1, much larger than has been seen in any radio pulsar. 

One possibility that can explain the large G for IE 2259+586, the large transient in- 
creases in the magnitude of z> in all three large glitches, as well as the large extended z> 
change in the first IE 1841—045 glitch, is that core superfluid is somehow involved, as it 
is expected to carry the bulk of the moment of inertia. We note that core glitches have 
been disc ussed in the radio p ulsar context for some time, albeit for very different reasons. 
Although lAlpar et al.l ( 1984bf ) argu ed that the crust and core should be strongly coupled on 
very short time scales, Jones! ( 1998 ) found that crustal pinning of superfluid vortices cannot 
be occuring in neutron stars, because the maxim um pinning force is orders o f magnitude 
smaller than the estimated vortex Magnus force. iDonati fc Pizzocherd (120031 ) argue that 
crustal vortex pinning cannot occur for independent reasons. If these authors are correct, 
pulsar glitches would generally not originate in the crust. 

Why would the clearest evidence for core glitches come from AXPs? The interaction 
between vortices and quantized magn etic flux tubes in a core superfluid could provide resis- 
tance to outward motion of vortices (jJonea Il998l ; iRuderman et al.1 11998c I Jonesl 120021 ) . The 
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interior magnetic field playing a role in vortex pinning as studied by iRuderman et all (119981 ) 
could help explain why such unusual glitch recoveries are seen preferentially in AXPs, which 
appear to have much larger magnetic fields than conventional radio pulsars. Perhaps the 
larger field, which implies a higher density of flux tubes, can effectively pin more superfluid 
vortex lines in a magnetar core; with greater magnetic activity, sudden magnetic rec o nfigu- 
rati ons would result in lar ge core vortex re configurations. W e , no ^ that iKaspi et al.1 (120001 ) 
and iDairOsso et al.1 (120031 ) argued that the IRuderman et al.1 (119981 ) model must be inappli- 



cable to AXPs as that m odel predicts no 



itche s for periods greater than ~ 0.7 s. However 



a more careful reading of IRuderman et al.1 (119981 ) reveals that this prediction does not apply 
for magnetar-strength magnetic fields. 



As pointed out by IKaspi et al.l (120001 ). the high temperatures of AXPs, as measured 
from their X-ray spectra, are at odds with the glitch observations. This is because in the 
crustal pinning models, the pinning force is highly temperature dependent, such that vortex 



lines can creep outward much more easily when the temperature is high (e.g., lAlpar et al. 



19891 ). This has long been the explanation (e.g., I Anderson fc Itohlll975l ) for the difference 
between the Crab and Vela pulsar glitch behaviors: the hotter Crab pulsar glitches less 
frequently and with smaller frequency jumps because its vortex array can move outward 
more smoothly. If this were true, AXPs, having measured effective temperatures much 
higher than the Crab pulsar, should glitch less frequently and with smaller glitch amplitudes 
than the Crab pulsar, clearly not what is observed. If we abandon the crustal glitch model 
(at least in AXP glitches) the absence of the expected temperature dependence a nd the 
observ ed universal age correlation could be explained. For example, as discussed by [Jones 
(119981 ). perhaps relatively smooth outward motion of core vortices could take place before 
the magnetic flux distribution necessary for impeding them has developed. In this picture, 
magnetic field distribution, not temperature, is the age- associated property that is a primary 
factor determining glitch behavior. 



Finally. ICampana et al.1 (120071 ) argued, on the basis of seven observations of RXS J170849. 0—400910 
obtained over ~10 yr, that observed spectral and flux variations were correlated with glitch 
epoch. They predict, given apparent flux increases seen in mid-2004 and mid-2005, that a 
glitch should occur soon thereafter (see Fig. [T4T) . Indeed as we have shown (see Table H]), 
an unambiguous glitch occured just following their mid-2005 observation. On the other 
hand, the first candidate glitch (Table [5]) occured when the total flux was very low and 
apparently declining. Thus, if there is a causal connection between long-term flux variations 
in RXS J170849. 0—400910 and glitches, either this candidate glitch is not a true glitch, 
or accompanying radiative changes are only relevant to large glitches. The sparsity of the 
total flux measurements, along with the relatively short time scale of the pulsed flux vari- 
ations we report in RXS J170849. 0—400910, suggest that more dense total flux monitoring 
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could reveal yet unseen fluctuations that are not glitch-associated. The flux variability of 
RXS J170849.0-400910 is discussed further below. 



6.2. Radiative Changes 



The approximate stabilities of the pulsed fluxes of RXS J170849.0-400910 and IE 1841-045 
fsee © are in contrast to those seen for AXPs IE 2259+586 and IE 1048.1-5937, both 



of which have shown large pulsed flux variations (IKaspi et all 120031 ; I Woods et all 12004 



Gavriil fc 



with time ( 



K^aspil I2004T) . and even 4U 0142+61 which has shown a slow pulsed flux increase 



Dib et al.ll2007ah . It is interest ing that the phase-averaged flux of RX S J170849.0-400910 



has been reported to be highly variable (IRea et al.ll2005l ; iCampana et al.1 120071 and Fig. [H]) . 
with changes as large as ~60% in 2004-2005, while the pulsed flux is not, with maximum 
contemporaneous change of < 15%. Note that this conclusion holds even when the phase- 
averaged fluxes, reported in the 0.5-10 keV band, are converted to 2-10 keV, that used 
for our RXTE observations. This suggests an anti-correlation between pulsed fraction and 
total flux that acts to ensure that the pulsed flux is roughly constant. If so, pulsed flux 
is not a good indicator of total energy output for RXS J170849.0— 400910A similar anti- 
correlation between total flux and pulsed fraction has been reported for AXP IE 1048—5937 



(ITiengo et al.ll2005l : I Gavriil et al 



flux does not remain constant (IGavriil fc Kaspi 



2006 



Gavriil et al.ll2007f ). although in that case the pulsed 



20041 ). but follows the phase-averaged flux, 



just with lower dynamic range (Tarn et al., submitted). For RXS J170849. 0—400910, the ex- 
actness of the anti-correlation is perhaps surprising. It could be that all the phase-averaged 
flux variations are in the 0.5-2 keV band, invisible to RXTE. However this woul d not jibe 



with t he reported correlatio n of the phase-averaged fluxes with power-law index (iRea et al. 



20051 ; ICampana et al.ll2007l ). We also note that in the phase-averaged flux analysis, the 
equivalent hydrogen column Nh was allowed to vary from observation to observation, rather 
than being held fixed at a constant. This inconsistency could bias the comparison, though 
likely not by a large amount. It is tempting to question the relative calibrations of the differ- 
ent instruments used to measure the phase-averaged flux of RXS J170849. 0—400910, as the 
greatest dynamic range is implied by lone XMM-Newton and Swift observations, and even 
the two Chandra X-ray Observatory observations were obtained with different instruments. 
Still, admittedly, relative systematic calibration uncertainties are not expected to yield a 
>50% dynamic range, as is reported. Regular monitoring with a single imaging instrument 
could settle this issue. 



Changes in pulse profile seem to be generic in AXPs and, as discussed above, are not 
always correlated with glitches, although in some cases, e.g., IE 2259+586, 4U 0142+61, 
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and possibly RXS J170849. 0—400910, they are. Spectrally no clear pattern ha s emerged. In 



IE 22 59+586, the pulse profil e changes followin g: its 2002 event were broadband (IWoods et al. 



2004 ). while in 4U 0142+61 flDib et al.ll2007al ) and IE 1841-045, they appear more promi- 
nent at soft energies. In RXS J170849. 0—400910, the changes are more apparent in the 
hard band. In the context of the magnetar model, this hints at crustal motions and surface 
activity, possibly coupled with magnetospheric activity, with the exact observational mani- 
festation dependent on a variety of factors ranging from viewing geometry to magnetospheric 
scattering optical depth. Whether ultimately this specific phenomenon will provide insights 
into the physics of magnetars remains to be seen. 



7. Summary 

We have reported on long-term RXTE monitoring of AXPs RXS J170849.0-400910 and 
IE 1841—045, which has allowed us to study these sources' timing, pulsed flux, and pulse 
profile evolutions. 

We have discovered four new AXP spin-up glitches, one in RXS J170849. 0—400910 and 
three in IE 1841-045, plus three new glitch candidates in RXS J170849.0-400910. Nearly 
all of the "classical" AXPs have now been seen to glitch, clearly demonstrating that this 
behavior is generic to the class. Moreover, in terms of fractional frequency increases, AXPs 
are among the most actively glitching neutron stars known. Further, unlike radio pulsar 
glitches, AXP glitches appear to come in two varieties: those that, like radio pulsars, are 
radiatively quiet in pulsed flux, and those that are, unlike radio pulsars, radiatively loud, 
including correlated sudden flux increases and pulse profile changes. Thus far there is no clear 
correlation between AXP glitch size and whether or not it will be radiatively loud or quiet - 
two of the largest AXP glitches thus far were quiet. We have found a substantial long-term 
increase in the magnitude of the spin-down rate in the largest glitch from IE 1841—045, and 
have also shown that large AXP glitches often have recoveries that are unusual compared 
with those seen in radio pulsars. Specifically, their spin-down rates in the days and sometimes 
weeks after a glitch are significantly larger in absolute value than their long-term spin-down 
rate. This latter effect may indicate a temporary reversal in the crust /crustal superfluid lag 
at the time of the glitch, or possibly more plausibly, and certainly more intriguingly, glitches 
of the core, which could explain the transient and extended v increases, as well as the large 
G value for IE 2259+586. 

Radiatively, we have found that the pulsed fluxes of RXS J170849.0-400910 and IE 1841-045 
are both fairly steady with time. This is perhaps surprising in light of the large changes in 
phase-averaged flux that have been reported for RXS J170849. 0—400910, and suggests, un- 
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less the latter are affected by systema tic calibration unce rtainties, that pulsed flux for this 



source, as for AXP IE 1048.1-5937 fjTiengo et al.l 120051 ). is not a good indicator of AXP 



X-ray output. Also, the pulse profiles of RXS J170849.0-400910 and IE 1841-045 both 
evolve; such evolution appears also to be a generic property of AXPs. However, no clear 
patterns in AXP pulse profile changes have yet emerged beyond occasional correlation with 
glitches. Hopefully further monitoring will shed physical light on this phenomenon. 
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Fig. 1— Epochs of observations of RXS J170849.0-400910 with RXTE. Gaps near the 
end/start of each year are due to Sun avoidance. See Tabled] for details. 
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Fig. 2. — Epochs of observations of IE 1841—045 with RXTE. Gaps near the end/start of 
each year are due to Sun avoidance. See Table [2] for details. 
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Fig. 3. — Spin and pulsed flux evolution in RXS J170849. 0—400910. Panels are described 
from top to bottom. Top: Frequency evolution, with inter-glitch intervals indicated for 
correspondence with ephemerides given in Table El Arrows indicate intervals for which glitch 
ephemerides were obtained (see Table Hj). Next: residuals, after subtraction of the best-fit 
models given in Table[3](with arbitrary inter-interval phase offsets subtracted). The increased 
scatter after MJD 52600 is due to a decrease in typical integration time and an increase in 
monitoring frequency. Next: Solid curve: frequency evolution of the models shown in Table [3] 
after removal of the linear trend defined by the frequency and frequency derivative from 
interval C as measured by fitting only those parameters. Data points: measured frequencies 
in independent sub-intervals after subtraction of the extrapolation of the same linear trend. 
Next: Evolution of the frequency derivative in sub-intervals, when fitting locally for only v 
and v. Bottom: Pulsed flux in the 2-10 keV range. All panels: Unambiguous glitch epochs 
are indicated with solid vertical lines. Candidate glitch epochs are indicated with dashed 
vertical lines. 
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Fig. 4. — The three unambiguous glitches observed in RXS J170849. 0—400910. Top panels: 
Frequency evolution around glitch as determined from ephemerides in Table SJ with the 
blow-up inset displaying glitch amplitude on a common scale for comparison. Middle panels: 
Residuals after subtraction of the best-fit pre-glitch ephemeris given in Table HI Bottom 
panels: Residuals after subtraction of best-fit glitch models given in Table HI All panels: 
Dashed vertical lines indicate assumed glitch epochs. 
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Fig. 5. — The three candidate glitches in RXS J170849. 0—400910. Top panels: Residuals 
after subtraction of best-fit pre-'glitch' ephemeris given in Table Middle panels: Residuals 
after subtraction of best-fit glitch models given in Table El Bottom panels: Residuals after 
subtraction of best-fit alternative models given in Table El All panels: Dashed vertical lines 
indicate assumed glitch epochs. 
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Fig. 6. — Spin and pulsed flux evolution in IE 1841—045. Panels are described from top to 
bottom. Top: Frequency evolution, with inter-glitch intervals indicated for correspondence 
with ephemerides given in Table [7J Arrows indicate intervals for which glitch ephemerides 
were obtained (see Tabled]). Next: Residuals, after subtraction of the best-fit models given in 
Tabled Next: Solid curve: frequency evolution of the models shown in Tabled after removal 
of the linear trend defined by the frequency and frequency derivative from the last year of 
data before the first glitch, as measured by fitting only those parameters. Dotted curve: 
alternate glitch recovery (see Section T3.2I for details). Filled circles: Measured frequencies 
in independent sub-intervals after subtraction of the extrapolation of the same linear trend. 
Unfilled squares: Epochs of the two immediate post-glitch observations (too few for the 
measurement of an independent frequency but crucial for the phase-coherent analysis). Next: 
Evolution of the frequency derivative in sub-intervals, when fitting locally for only v and v. 
Bottom: Pulsed flux in the 2-10 keV range. All panels: Glitch epochs are indicated with 
solid vertical lines. The dashed vertical line indicates the start of ephemeris B2, which does 
not include the two immediate post-glitch observations (indicated with unfilled squares). 
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Fig. 7. — The three glitches in IE 1841—045. Top panels: Frequency evolution around 
glitches as determined from ephemerides in Table [BJ with the blow-up inset displaying glitch 
amplitude. Middle panels: Residuals after subtraction of the best-fit pre-glitch ephemeris 
given in Table [HI Bottom panels: Residuals after subtraction of best-fit glitch models given 
in Table El All panels: Dashed vertical lines indicate assumed glitch epochs. 



31 



0.8 
0.4 


0.8 

0.-1 


O.ti 
0.4 



m 

a o 

0) 

B, 0.8 

s. 







05 



o 

Z 04 



0.8 
0.-1 


0.8 
0.4 




1 L II 1 1 II 1 LI 1 1 II 

J\ ' A ' (A) 

-i i i i 1 i r^t 1 1 1 i i 1 i vf\ 


A 
\ 


A 




(A) 






V w 




(CI) 

V V 






(CI) 












jA ' " " jA ' '(C2) 

\Aw 










' ' ' 


(D) 










(FG, 

- V v 


i ji i I i i i i I i ji i I i i i 





0.4 



0.8 
0.4 



0.4 



0.4 



OH 
0.4 



0.8 
0.4 



0.8 
0.4 



0.5 1 1.5 20 0.5 1 1.5 20 0.5 1 1.5 2 

Pulse Phase (cycles) Pulse Phase (cycles) Pulse Phase (cycles) 
2-10 keV 2-4 keV 4-10 keV 



Fig. 8. — Normalized pulse profiles in three energy bands for RXS J170849. 0—400910 for 
the seven glitch-free intervals (with corresponding labels at the top right) defined in the top 
panel of Figure [31 Different data qualities within an energy range are due to different net 
exposure times. Two cycles are shown for each profile for clarity. 
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Fig. 9. — Left panel: Time evolution of the ratio of the power in the nth harmonic to the 
total power in the 2-10 keV pulse profile of RXS J170849.0-400910. Circles represent n = 1, 
squares n — 2, and triangles n — 3. Solid vertical lines indicate epochs of glitches; dashed 
vertical lines are epochs of candidate glitches. Middle panel: Same as left panel but for 
2-4 keV. Right panel: Same as middle but for 4-10 keV. 
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Fig. 10. — Top panel: Time evolution of the ratio of the power in the nth harmonic to the 
total power in the 2-10 keV pulse profile of RXS J170849. 0—400910. The circles represent 
n — 1, squares n — 2, and triangles n — 3. Solid and dashed vertical lines indicate epochs of 
glitches and candidate glitches, respectively. The probability that the observed fluctuations 
are due to random noise are 68%, 97% and 69% for n — 1, 2, 3, respectively. Bottom panel: 
reduced \ 2 P er degree of freedom for successive profile differences (see text §4. II for details). 
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Fig. 11. — Normalized pulse profiles in three energy bands for IE 1841—045 for the five 
glitch-free intervals defined in the top panel of Figure Different data qualities in each 
energy range are due to different net exposure times. Two cycles are shown for each profile 
for clarity. 
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Fig. 12. — Left panel: Time evolution of the ratio of the power in the nth harmonic to the 
total power in the 2-10 keV pulse profile of IE 1841—045. Circles represent n — 1, squares 
n — 2, and triangles n — 3. Solid vertical lines indicate epochs of glitches. Middle panel: 
Same as left panel but for 2-4 keV. Right panel: Same as middle but for 4-10 keV. 
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Fig. 13. — Top panel: Time evolution of the ratio of the power in the nth harmonic to 
the total power in the 2-10 keV pulse profile of IE 1841—045. Circles represent n — 1, 
squares n = 2, and triangles n = 3. Solid vertical lines indicate epochs of glitches. The 
probabilities that the observed fluctuations arise from random noise are 99%, 97% and 96% 
for n = 1, 2, 3, respectively. Bottom panel: reduced \ 2 P er degree of freedom for successive 
profile differences (see text §4.21 for details). 
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Fig. 14. — Frequency, pulsed flux, reported total unabsorbed flux, and reported photon 
index as a function of time for RXS J170849. 0—400910. Frequency and pulsed flux data 
are identical to those shown in Figure [31 Solid and dashed vertical lines indicate epochs of 
glitches and glitch candidate s, respectively. U nabs orbed phase-ave r aged 0.5-10 keV fluxes 
and photon indexes are from iRea et al.l (120051 ) and ICampana et al.l ( 120071 ) , and are labelled 
by observing telescope. That the pulsed flux remains relatively constant while the phase- 
averaged flux appears to vary by nearly a factor of two (albeit as measured by different 
instruments) suggests a strong anti-correlation between total flux and pulsed fraction. 



10 -8 -6 -4 

LOG[Az//V] 




-10 -8 -6 -4 

LOG[Az/] 



Fig. 15. — Amplitude distribution of AXP glitches (bold line) and radio pulsar glitches (thin 
line) for (a) fractional frequency jump and (b) absolute frequency jump (in Hz). Radio pulsar 
glitch amplitudes are from an unpublished catalog kindly supplied by A. Lyne. AXP glitches 



included here are those listed in Tables 0] and El the 2002 IE 2259+586 glitch (IKaspi et al. 



20031 : IWoods et al.ll2004f ). as well as a recent unpublished IE 2259+586 glitch (Dib et al., in 
prep.). 
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Fig. 16. — Activity parameters versus age (as estimated from vjlv) and versus z> for radio 
pulsars and AXPs. Fractional activity a g is defined using the sum of the fractional frequency 
changes, while activity A g is defined using the absolute frequency jumps. The only radio 
pulsars included (filled circles) are those having exhibited three glitches or more during 
continual (e.g., bi-monthly) monitoring, as recorded in the unpublished glitch catalog kindly 
supplied by A. Lyne. The AXPs included here (open circles) are RXS J17084 9.0— 400910, 



IE 1841-045, and IE 2 259+586. The latter has glitched twice, once in 2002 (IKaspi et al. 



20031 ; IWoods et al.ll2004l ). and once in 2007 (Dib et al., in prep.). Only unambiguous glitches 



were included for RXS J170849. 0—400910; as the candidate glitches are small, including 
them does not make a qualitative difference. 
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Fig. 17. — "Coupling parameter" G (as defined by (ILink et al.lll999l ); see §6.11) as a function 
of spin-down age (z//2z>). Solid points are radio pulsars with three or more observed glitches, 
as recorded in the unpublished catalog of A. Lyne. Open circles are AXPs IE 2259+586, 
RXS J170849.0-400910, and IE 1841-045. 
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Table I. Summary of RXTE Observations of RXS J170849.0-400910 



Obs. 


Typical 


Typical 


No. of 


Total 


First - Last 






Cycle 


Exp. a 


Separation 3, 


Obs. b 


Exp. c 


MJD d 


First Date - 


Last Date 




(ks) 


(days) 




(ks) 








3 


2.5 


15 


29 


75 


50825.7-51186.7 


01/12/1998- 


-01/08/1999 


4 


3 


24 


20 


60 


51215.7-51614.1 


02/06/1999- 


-03/11/2000 


5 


3 


24 


13 


40 


51655.7-52041.5 


04/21/2000- 


-05/12/2001 


6 


3 


29 


13 


40 


52049.5-52325.6 


05/20/2001- 


-02/20/2002 


7 


5.5 


23 


12 


65 


52366.5-52718.7 


04/02/2002- 


-03/20/2003 


8 


1.8 


5 


58 


105 


52745.7-53058.6 


04/16/2003- 


-02/23/2004 


9 


2 


5 


70 


135 


53063.1-53429.1 


02/28/2004- 


-02/28/2005 


10 


2 


8 


47 


90 


53435.1-53791.6 


03/06/2005- 


-02/25/2006 


ll c 


2 


7 


32 


60 


53799.0-54015.4 


03/05/2006- 


-10/07/2006 



a The exposure and separation are approximate. Note that the PCA effective area changed with time 
primarily due the reduction of the average number of PCUs operational during an integration. This effect 
is not incorporated in the tabulated integration times. 

b When the last digits of the observation ID of two successive data sets are different, the two data sets are 
considered separate observations. 

c The total exposure does not include Earth occultation periods. 

d First MJD and Last MJD are the epochs, in Modified Julian Days, of the first and the last observations 
in a Cycle. 

c Cycle 11 not yet completed. 
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Table 2. Summary of RXTE Observations of IE 1841-045 



Obs. 
Cycle 


Typical 
Exp. a 
(ks) 


Typical 
Separation 3 - 
(days) 


No. of 
Obs. b 


Total 
Exp. c 
(ks) 


First - Last 
MJD d 


First Date - 


- Last Date 


4 


4.5 


27 


26 


120 


51224.4-51597.3 


02/15/1999- 


-02/23/2000 


5 


4.5 


38 


16 


70 


51644.7-51976.9 


04/10/2000- 


-03/08/2001 


6 


7 


27 


8 


50 


52001.6-52300.1 


04/02/2001- 


-01/26/2002 


7 


12 


45 


7 


80 


52349.8-52666.0 


03/16/2002- 


-01/27/2003 


8 


4.5 


20 


17 


80 


52726.8-53052.9 


03/28/2003- 


-02/17/2004 


9 


4.5 


20 


19 


80 


53073.7-53413.3 


03/09/2004- 


-02/12/2005 


10 


5 


14 


31 


130 


53440.0-54153.9 


03/11/2005- 


-02/22/2007 


ll c 


5 


14 


12 


60 


53800.9-53970.6 


03/06/2006- 


-08/23/2006 



a The exposure and separation are approximate. Note that the PCA effective area changed with time primarily 
due the reduction of the average number of PCUs operational during an integration. This effect is not incorporated 
in the tabulated integration times. 

b When the last digits of the observation ID of two successive data sets are different, the two data sets are 
considered separate observations. 

c The total exposure does not include Earth occultation periods. 

d First MJD and Last MJD are the epochs, in Modified Julian Days, of the first and the last observations in a 
Cycle. 

c Cycle 11 not yet completed. 
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Table 3. Long-Term Spin Parameters for RXS J170849.0-400910 a 



Parameter 




Ephcmcris A 
Spanning MJD 
50826-51418 


Ephcmcris B 
Spanning MJD 
51446-51996 


Ephcmcris C 
Spanning MJD 
52036-52960 


Ephcmcris D 
Spanning MJD 
53010-53325 


Ephcmcris E 
Spanning MJD 
53377-53548 


Ephcmcris F 
Spanning MJD 
53556-53631 


Ephcmcris G 
Spanning MJD 
53638-54015 


MJD start 




50826.078 


51446.610 


52035.655 


53010.094 


53377.133 


53555.734 


53638.033 


MJD end 




51418.374 


51995.680 


52960.186 


53325.061 


53547.811 


53631.161 


54015.487 


TOAs 




39 


19 


74 


69 


29 


13 


49 


v (Hz) 




0.090913818(2) 


0.090906071(3) 


0.090906089(3) 


0.090892731(13) 


0.090887608(18) 


0.090885281(8) 


0.090884082(9) 


v (10~ 13 Hz s -1 ) 




-1.583(3) 


-1.574(2) 


-1.565(6) 


-1.40(5) 


-1.19(8) 


-1.70(2) 


-1.58(3) 


v (10~ 22 Hz s~ 2 ) 




-1.4(3) 


0.36(9) 


-8.9(1.8) 


-44(11) 


-131(23) 




-8(7) 


d 3 v/dt 3 (10~ 28 Hz s 


" 3 ) 


-0.056(9) 




1.5(3) 


5.5(1.6) 


15(3) 




1.4(9) 


d 4 i//dt 4 (10~ 35 Hz s 


- 4 ) 






-1.4(3) 


-3.1(1.0) 






-0.8(6) 


d 5 v/dt 5 (10~ 43 Hz s 


" 5 ) 






7.4(1.5) 










d 6 v/dt 6 (10~ 5Q Hz s 


-6j 






-1.8(4) 










Av d b (Hz) 








36(3) X10~ OS 










t d b (days) 








43(2) 










Epoch (MJD) 




51445.3846 


52016.48413 


52016.48413 


52989.8475 


53366.3150 


53549.15095 


53635.6772 


RMS residual (phase) 




0.0079 


0.0150 


0.0154 


0.0132 


0.0142 


0.0112 


0.0154 



a Numbcrs in parentheses are TEMPO-rcportcd lcr uncertainties. 

^Parameters held fixed at values determined from local glitch fits as shown in Tableful 
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Table 4. Local Ephemerides of RXS J170849.0-400910 Near Glitch Epochs 3 





Ephemeris 


Ephemeris 


Ephemeris 


Parameter 


Near Glitch 1 


Near Glitch 2 


Near Glitch 3 



MJD range 
TOAs 

Epoch (MJD) 
v (Hz) 

v (1(T 13 Hz s- 1 ) 
Glitch Epoch (MJD) 
Ai/ (Hz) 
Az> (Hz s- 1 ) 
Au d (Hz) 
t d (days) 

RMS residual (phase) 



51186.503-51614.187 51614.185-52366.663 53465.392-53631.161 



22 

51445.3846 
0.090913822(2) 
-1.5714(14) 
51445.3846 
5.1(3)xl0~ 8 
-0.8(4)xl0~ 15 



0.0102 



29 

52016.48413 
0.090906068(2) 
-1.5797(11) 
52016.48413 
2.2(4) xl0~ 8 
-l.l(2)xl0- 15 
36(3) xl0~ 8 
43(2) 
0.0193 



26 

53549.15095 
0.090885035(9) 
-1.67(2) 
53549.15095 
24.6(9) xl0~ 8 
-2(2)xl0~ 15 



0.0140 



''Numbers in parentheses are TEMPO-reported la uncertainties. 
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Table 5. Local Ephemerides of RXS J170849.0-400910 Near Candidate Glitch 

Epochs a 



Parameter 


Ephemeris 
Near Candidate 1 


Ephemeris 
Near Candidate 2 


Ephemeris 
Near Candidate 3 


MJD range 


52745.790-53140.604 


53229.271-53456.688 


53562.209-53785.6E 


TOAs 


78 


38 


28 


Epoch (MJD) 


529898475 


53366.3150 


53635.6772 


v (Hz) 


0.0908927493(18) 


0.090887617(5) 


0.090884020(8) 


v (1(T 13 Hz s- 1 ) 


-1.5842(15) 


-1.570(7) 


-1.67(2) 


Glitch Epoch (MJD) 


529898475 


53366.3150 


53635.6772 


Az> (Hz) 


2.8(4)xl0~ 8 


5.2(6) xl0~ 8 


6.7(3) xl0~ 8 


Az> (Hz s- 1 ) 


b 


-1.9(1.3)xl0- 15 


6.0(5)xl0- 15 


Au d (Hz) 








t d (days) 








RMS residual (phase) 


0.0153 


0.0099 


0.0110 



Numbers in parentheses are TEMPO-reported la uncertainties. 
'Entries with the value '0' are consistent with being zero. 
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Table 6. Alternate Ephemerides of RXS J170849.0-400910 Near Candidate 

Glitch Epochs a 



Parameter 


Ephemeris 
Near Candidate 1 


Ephemeris 
Near Candidate 2 


Ephemeris 
Near Candidate 3 


MJD range 


52745.790-53140.604 


53229.271-53456.688 


53562.209-53785.6E 


TOAs 


75 


38 


28 


Epoch (MJD) 


52989.8475 


53366.3150 


53635.6772 


v (Hz) 


0.0908928173(6) 


0.0908876402(13) 


0.090884059(2) 


z> (10~ 13 Hz s" 1 ) 


-1.556(3) 


-1.492(6) 


-1.485(16) 


v (1(T 22 Hz s" 2 ) 


1.2(4) 


4(2) 


b 


d 3 is/dt 3 (10~ 28 Hz s~ 3 ) 


-0.66(16) 


-6.6(1.4) 


-19(4) 


d 4 u/dt 4 (lO" 35 Hz s~ 4 ) 


-1.0(4) 


-18(6) 


52(12) 


RMS residual (phase) 


0.0146 


0.0106 


0.0152 



Numbers in parentheses are TEMPO-reported la uncertainties. 
'Entries with the value '0' are consistent with being zero. 
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Table 7. Long-Term Spin Parameters for IE 1841— 045 a 



Parameter 




Ephcmcris A 
Spanning MJD 
51225-52438 


Ephcmcris B 
Spanning MJD 
52460-52981 


Ephcmcris B2 
Spanning MJD 
52610-52981 


Ephcmcris C 
Spanning MJD 
53030-53816 


Ephcmcris D 
Spanning MJD 
53829-53983 


MJD start 




51224.538 


52460.000 


52610.313 


53030.093 


53828.808 


MJD end 




52437.712 


52981.186 


52981.186 


53815.842 


53983.431 


TOAs 




53 


19 


17 


54 


11 


v (Hz) 




0.0849253002(9) 


0.084904428(7) 


0.084889922(3) 


0.084890135(6) 


0.084868767(7) 


v (10~ 13 Hz s -1 ) 




-2.9940(10) 


-3.176(2) 


-3.179(2) 


-3.354(7) 


-2.833(9) 


v (10~ 22 Hz s~ 2 ) 




3.30(14) 






16.4(4) 




d 3 u/dt 2 (10~ 29 Hz s 


" 3 ) 


0.9(2) 






-2.81(11) 




d i v/dt i (10 -36 Hz s' 


-") 


-2.3(2) 










d 5 u/dt 5 (10~ 43 Hz s 


" 5 ) 


1.5(3) 










d 6 v/dt 6 (10~ 51 Hz s 


-B) 


-3(2) 










Av d b (Hz) 






8.1(6)xl0~ 7 








t d b (days) 






43(3) 








Epoch (MJD) 




51618.000 


52464.00448 


52997.0492 


52997.0492 


53823.9694 


RMS residual (phase) 




0.029 


0.025 


0.028 


0.033 


0.022 



a Numbcrs in parentheses arc TEMPO-rcportcd lcr uncertainties. 

^Parameters held fixed at values determined from local glitch fits as shown in Table [Fl 
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Table 8. Local Ephemerides of IE 1841-045 Near Glitch Epochs 3 





Ephemeris 


Ephemeris 


Ephemeris 


Parameter 


Near Glitch 1 


Near Glitch 2 


Near Glitch 3 



MJD range 
TOAs 

Epoch (MJD) 
v (Hz) 

v (1(T 13 Hz s- 1 ) 
Glitch Epoch (MJD) 
Ai/ (Hz) 
Az> (Hz s- 1 ) 
^v d (Hz) 
t d (days) 

RMS residual (phase) 



52001.684-52981.186 52773.845-53244.179 53579.360-53970.681 



30 

52464.00448 
0.084903950(2) 
-2.8980(10) 
52464.00448 
4.78(7) xl0~ 7 
-2.78(2)xl0" 14 
8.1(6)xl0~ 7 
43(3) 
0.022 



26 

52997.0492 
0.084889815(3) 
-3.162(3) 

52997.0492 
2.08(4) xl0~ 7 

4(3) xl0~ 16 



0.015 



31 

53823.9694 
0.084868657(4) 
-2.872(4) 
53823.9694 
1.18(7)xl0~ 7 
2(l)xl0~ 15 



0.022 



''Numbers in parentheses are TEMPO-reported la uncertainties. 



